Incompressible Quantum Glass State of Bosons in Two Dimensions 



Yanchcng Wang,^ Wcnan Guo,^'El and Anders W. Sandvik^'0 

'^Department of Physics, Beijing Normal University, Beijing 100875, China 
^Department of Physics, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, USA 

(Dated: October 17, 2011) 

We study the quantum glass state intervening between the conventional superfluid and Mott- 
insulator states of lattice bosons with random potentials. Its properties at temperature T = Q are 
controlled by rare large regions of superfluid surrounded by Mott insulator. These regions make the 
state gapless although it is insulating. Contrary to the commonly accepted theory of this state in two 
dimensions, we show here that a vanishing gap does not necessarily imply nonzero compressibility. 
Using quantum Monte Carlo simulations of the Bose-Hubbard model and a percolation theory, we 
show that the compressibility k at average filling p = 1 follows the form k ~ exp(— 6/T") with 
a < 1. The system is, thus, incompressible at T = and should be classified as a Mott glass. 

PACS numbers: 67.85.-d, 67.85.Hj, 67.10.F, 64.70.Tg 



Interacting lattice bosons have two distinct types of 
ground states in the absence of disorder; the superfluid 
(SF) and the Mott-insulator (MI). In the standard proto- 
typical Bose-Hubbard (BH) model with repulsive on-site 
interactions, an MI state has an integer number of par- 
ticles per site and there is a gap to states with added or 
removed particles. The gapless SF can have any filling 
fraction. These phases and the quantumphase transi- 
tions between them are well understood [iHa] and have 
also been realized experimentally with ultracold atoms 
in optical lattices @, 01 • If disorder in the form of ran- 
dom site potentials is introduced in the BH model a third 
kind of state appears — an insulating but gapless quantum 
glass (QG). This state has been the subject of numerous 
studies over the past two decades [J-^S, 8 20] , but many 
of its properties are still not well understood. In one di- 
mension, two typses of QG ground states are known; the 
compressible Bose glass (BG) and the incompressible 
Mott glass (MG) [l9(, with the latter appearing at com- 
mensurate filling fractions only [20| . In two and higher 
dimensions it has been believed that the QG is always 
of the compressible BG type [2, IJ. QG states are also 
accessible in optical lattices [2l|, [2^ ■ 

We here discuss site-disordered bosons at temperature 
T > 0. Our arguments and conclusions will mostly be 
generic, but for definiteness, to study systems using unbi- 
ased quantum Monte Carlo (QMC) simulations, we con- 
sider the two-dimensional (2D) BH Hamiltonian 
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where (ij) denotes nearest-neighbor sites on the square 
lattice, 6^ (6^) are boson creation (destruction) opera- 
tors, rii = bfbi site occupation numbers, and random 
potentials uniformly distributed in the range [— A — /i. A — 
fj] around the average chemical potential /i. 

Previous studies of the site-disordered BH model have 
focused on its T = phase diagram and, in particular, 



the long-standing question of whether or not a QG state 
always intervenes between the SF and MI phases in the 
presence of disorder 0, [sl-fisj . Such questions arise in 
many classical and quantum many-body problems, with 
the QG corresponding generically to a Griffiths phase 
in which statistically rare large regions of some phase A 
inside another phase B lead to singularities character- 
izing a new type of disordered phase (Griffiths phase) 
not present in the absence of disorder j23l - f25j . In the 
case of the BH model with bounded disorder (which is 
necessary to retain an MI state), the Griffiths argument 
states the following 0, Q : In the MI, once the width 2 A 
of the disorder distribution exceeds the Mott gap Am, 
there can be arbitrarily large regions of SF inside the 
MI. These SF "puddles" are not immediately, until A ex- 
ceeds some larger critical value, percolating through the 
lattice, however, so the state is insulating. The QG is 
still fundamentally different from the MI because it is 
gapless. This scenario was recently refined and named 
the "theorem of inclusion" [iJl ■ 

While, according to the theorem, the MI-QG bound- 
ary is dictated exactly by the Mott gap of the clean MI 
B 13 ' it tis-s been pointed out that this boundary can- 
not be extracted based on QMG simulations, because the 
rare regions leading to the vanishing gap are in practice 
not generated on the relatively small lattices to which 
the calculations are limited (1, 14, Indeed, some 

earlier numerical studies d, [lol | had partially missed the 
QG phase because its signatures at T = are difficult to 
distinguish from those of the MI. Renormalization group 
(RG) studies as well had pointed to the possibility of a 
direct MI-SF transition [lll[l3- 

The purpose of our work presented here is to clarify the 
consequences of gaplessness of the 2D QG — in particular, 
to show that the absence of a gap does not necessarily im- 
ply a compressible state. We will argue that, in contrast 
to the commonl y a cce p ted theory of 2D site-disordered 
bosons [liSEjIlS, [i^j the compressibility k of the QG 
can vanish when T 0. For the model llj at filling 
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FIG. 1: (Color online) Disorder-averaged compressibility ver- 
sus inverse temperature for systems with p = 1, U = 22, and 
different values of the width A of the disorder distribution. 
The lattice size is L = 32 in all cases. The curves are fits to 
the exponential form ([2]) with a = 1 for A = 0, 3 (MI state) 
and a = 0.77 and 0.53 for A = 6 and 7, respectively (MG 
state). For A = 9 (SF state) the line is horizontal. 

fraction p = 1 wc find the form 

n^exp{-b/T^), (2) 

where a < 1. In the MI phase the presence of a gap A 
and a slowly varying density of states at the gap edge 
implies the above form with a = 1 and b = A. The 
incompressibility of the QG at T = puts it in the MG 
class, which in site-disordered systems was up until now 
only found in one dimension [19[. 

We next present QMC results demonstrating ^ and 
then discuss some aspects of the density of states and 
how the anomalous exponential compressibility follows 
from an SF quantum percolation scenario. 

QMC calculations — We have used the stochastic series 
expansion QMC method with directed-loop updates [26 1 
to study the Hamiltonian ([1]) with the chemical potential 
fj. adjusted to give an average filling fraction p = (n) /N = 
1 (typically to accuracy 10~^ or better) in averages over 
hundreds of realizations of the random potentials. To 
speed up the simulations we restrict the site occupation 
to rii < 2, which for large U only has minor unimportant 
consequences (which we have checked explicitly by some 
calculations with rii < 3). We fix C//t = 22. According to 
previous work [13] (with which our results agree well) the 
system is then an MI for A < 4.2, an SF for A > 7.8, and 
a QG in the window between these points (and also in a 
separate region for A > 70 that we do not study here). 
For inverse temperatures /? = i/T up to 8, we eliminate 
finite-size effects by using Lx L lattices with L up to 32. 
Results for L = 16 and L = 32 agree within error bars. 

As shown in Fig. [U the behavior in the QG (MG) 
phase follows the form ([2|) with a « 0.77 for A = 6 and 
a ~ 0.53 for A = 7. For points inside the MI phase 
(A = and 3) the behavior agrees with the conventional 
a = 1 exponential as expected (and also discussed further 
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FIG. 2: (Golor online) Compressibility distribution for L = 16 
systems at different inverse temperatures. The total number 
of instances of the random potentials was 1200 for each /3. 

below). Inside the SF phase, exemplified in the figure by 
results for A = 9, the compressibility converges rapidly 
to a non-zero value as T is lowered, as also expected. 

An important question is whether the low-T k distri- 
bution develops a stretched tail, which would make it dif- 
ficult to obtain reliable disorder averages based on a rel- 
atively small number of samples. This problem at T = 
was stressed in 

Fig.Hshowshisto grams for L = 16 
systems. It is clear that the distribution becomes increas- 
ingly skewed as the temperature is lowered, with a thin 
high-K tail developing. The peak moves rapidly toward 
lower values, however, and the tail is not sufficiently long 
or fat to cause serious problems with properly converging 
the average k. We have confirmed that successive bin- 
ning of the samples into increasingly large groups leads 
to convergent error bars (which in Fig. [1] are smaller than 
the plot symbols in most cases). It is likely that the tails 
will become very significant at even lower temperatures, 
but with the whole distribution shifting toward lower val- 
ues, so that the average k — >■ as T — > according to 
the anomalous exponential form ([2|). 

Compressibility and density of states — According to 
Eq. ^ the QG is incompressible at T = 0, i.e., it is an 
MG, in contrast to the long-held belief that it should be 
a compressible BG 0, H, [H , 14 1 . The gaplessness in the 
thermodynamic limit follows from the presence of rare 
large regions of SF surrounded by MI. An infinite SF do- 
main is of course compressible and leads to a vanishing 
gap of the system as a whole. The probability of a given 
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site belonging to an infinitely large SF region is vanish- 
ingly small, however, and the compressibility normalized 
by the system volume in the standard way can therefore 
vanish when T — >■ despite the absence of a gap. 

We will next consider the basic relationships between 
the density of states and the compressibility, to demon- 
strate more formally why an MG with k — > is plausible. 
To this end, we write the compressibility as 

" TN E„/*e-/^p„(e) ' 

where Pni^) is the density of states with n particles. 

Consider first the MI state of a clean system. There 
is a gap A to adding a particle. In a simplified picture, 
valid for large U, the added particle should disperse as a 
free boson, i.e., with the 2D energy- momentum relation 
Eq = A + 4:t — 2t[cos{qx) + cos(gy)]. This gives an almost 
flat density of states at the band edges and a van Hove 
singularity at the center. At low temperatures the singu- 
larity should not matter and the compressibility should 
be governed by the density of states at the lower band 
edge. We can therefore consider a simplified flat density 
of states, which gives the same low-T behavior (up to 
unimportant T and N independent factors). 

For definiteness, we again consider a system with av- 
erage density p = 1. There are N single-particle states, 
for both added and removed particles, and the energies 
in these sectors are the same due to the emergent Q 
particle-hole symmetry. We distribute the states uni- 
formly within a band of width W above the gap A. 
Thus, the density of these states is Pn±i{^) = N/W for 
e € [A, A -I- W] and otherwise. The compressibility at 
low T, where we only need to consider n ^ N in the 
denominator of ©, is then k = (2/W^)c"'^/^. The pure 
exponential form agrees with our A = and A = 3 QMC 
results in Fig. [TJ For A > A should correspond to an 
average gap, smaller than the gap of the clean system. 

Note that, in a clean system with perfect particle-hole 
symmetry, the thermodynamic gap A controlling k as 
above should be equal to half the Mott gap Am- The 
A = fit in Fig. [Hgives A = 4.3(1), which indeed agrees 
well with Am w 8.34 extracted from the difference in 
chemical potential between the upper and lower edges of 
the p = 1 Mott lobe at J7 = 22. 

When the system is gapless it is no longer true that 
the ground state dominates the partition function at low 
temperatures. Consider the n — N zL 1 sectors. When 
A'^ —> cx) we can neglect the ground state contribution to 
the partition function in ([3]) if pN±i{i ^0) > or if 
P7V±i(f) vanishes sufficiently slowly when e — >■ 0. Since 
{n — N)'^ = 1 we get a compressibility contribution cx 
1 /NT, independently of the exact form of the density of 
states. Thus, the contribution from the single-particle 
sector vanishes for any T > as A^ — > oo. 

In the gapless case we can of course not neglect the sec- 
tors with more than one added or removed particle, be- 



cause those sectors should also be gapless, at least when 
the deviation from the ground state occupation number 
6n ^ N. In order for the total compressibility to be 
size independent (non-zero when A'^ — >■ oo) in the limit 
T — >■ 0, it is crucial, however, that there are enough low- 
energy states also in sectors with S„ up to order 0{y/N). 
This is required in order to cancel the factor 1/A^ in the 
normalization of k in ([3]). 

In the standard picture of the QG, a significant den- 
sity of low-energy states with 0{\/N) fluctuations seems 
unlikely, because the probability of a given site belong- 
ing to an SF region (which is where the particle number 
can fluctuate) decreases exponentially with the size of 
the region. Finite SF domains should have gaps between 
states in different particle number sectors, and, thus, be- 
come incompressible as T — >■ 0. These SF domain gaps 
will lead to a very low density of low-energy states for 
the sectors with n deviating significantly from N. The 
scenario k — > as T — > 0, suggested by our numerical 
results, therefore appears very plausible. 

Percolation scenario — The form Q can be understood 
heuristically using percolation theory: Consider an en- 
semble of SF domains below the percolation point. At 
T — l/m° all domains of sizes s < m are effectively in 
their ground states and do not contribute to the com- 
pressibility. The exponent a depends on the low-energy 
level spectrum of the SF domains, which may be affected 
by the domain shapes and boundaries. Domains with 
s > m should contribute essentially independently of T 
and s. The probability of a given site belonging to a 
domain with s > to is oc exp(— foTO'^), where c is known 
in classical percolation theory but may be different 
in a quantum system. In terms of the unknown expo- 
nents a and c the compressibility due to SF domains is 
K oc exp(-&T"'=/°), which is Eq. ^ with a = c/a. 

One might worry about the validity of this picture 
when the SF domains are not completely independent, 
being weakly coupled through the surrounding MI. How- 
ever, since the SF is not percolating through the system, 
the level structure should not be affected beyond a broad- 
ening of the individual domain levels on a much smaller 
scale than the level spacing. It remains to be understood 
how a can depend on the BH model parameters. 

Cross-over behavior — We have not yet carried out a 
systematic study of the BH model on a dense grid of dis- 
order strengths A, but from the results at two different 
MG points in Fig. [T] we know that the exponent a de- 
pends on the system parameters. There should always 
be a contribution oc exp(— A/T) to the compressibility 
from the MI background (with A an average gap), but 
this becomes negligible at low T when a < 1. Upon 
approaching the MG-MI transition, 6 in ([2]) should di- 
verge. The relative MI contribution then grows, so that 
K crosses over to the form with a = 1. At the T = 
MG-SF transition one expects k ^ {Ac - A)''^^""') with 
z = 2 (as in the generic clean MI-SF transition) 0. As 
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the transition is approached from the MG, a — ^ would 
ensure that k can evolve smoothly for any T > between 
the MG form and the constant k exactly at the transi- 
tion (with a discontinuity developing as T — > 0). Both 
the transitions out of the MG should be studied in more 
detail to clarify the exact nature of the cross-overs. The 
behavior at the MG~SF transition, in particular, should 
shed light on the mechanism of quantum percolation of 
the SF puddles. 

In a previous QMG study of a (2-1- l)-dimensional 
classical model with a glassy phase, mimicking the 2D 
quantum BH model ([T]), it was claimed that the T = 
compressibility at integer filling is non-singular at the SF 
transition [l3l | and nonzero in the glass phase. However, 
if a —> at this transition (and also in other possible 
cross-over scenarios), simulations carried out at low but 
non-zero temperature in the MG in close vicinity of the 
SF boundary may produce results that appear to con- 
verge to a non-zero value, even though k — >■ eventually 
as r — >■ according to ([2]) with a small a. In [l^j k > 
was indeed seen only in simulations very close to the SF. 

Conclusions and discussion — We have presented plau- 
sible arguments and unbiased QMG results demonstrat- 
ing an incompressible MG phase in the BH model. Previ- 
ous arguments for a generically compressible BG, starting 
with the seminal work by Fisher et al. Q and continu- 
ing up to more recent studies 3, have never, to our 
knowledge, been backed up by complete proofs taking 
into account the probability of 0{^/N) particle-number 
fluctuations. The outlines of proofs presented in Refs. [3 
and seem to neglect the fact that the normalization 
by in Eq. Q makes k — )- possible even though 
there are SF domains with arbitrarily small gaps. 

Our key physical argument is that a finite SF domain 
surrounded by an MI in the BH model becomes incom- 
pressible as T 0. This is clear for a single SF domain 
in a large MI, where the ground state has a definite num- 
ber of particles and there must be a gap to states with 
added or removed particles. The gap decreases with in- 
creasing size of the SF domain. In the QG state, with 
many SF domains surrounded by MI, as long as the SF 
cannot percolate through the system (which would bring 
it into the SF state) it should effectively consist of iso- 
lated SF puddles with fixed particle numbers at T — 0. 
If this picture is correct the anomalous exponential form 
([2]) follows, independently of dimensionality. 

While we have carried out substantial QMG calcula- 
tions only at density p ~ 1 (in part because it is much 
more difficult to locate QG states at significantly lower 
fixed average density), there is no apparent reason, based 
on the percolation scenario, why QG states of the BH 
model with non-integer average density should be quali- 
tatively different. In principle it is possible, however, that 
the emergent particle-hole symmetry Q at integer p in 
the clean system survives in the presence of disorder (as 
it can do in one dimension 



tems special. Indeed, according to a recent weak-disorder 
RG calculation [l^, 2D systems with integer filling have 
a different fixed point that is not accessible within that 
RG (and nothing could be said about the compressibil- 
ity) . Previous works 0, Hi found no distinction between 
integer and noninteger filling. Our results at p = 1 may 
then correspond to this unknwon fixed point, and the 
possibility remains of a compressible BG existing at non- 
integer filling. This would also imply that the percolation 
picture does not apply at noninteger fillings. 

The percolation scenario can, in fact, break down in 
one dimension, where, according to strong-disorder RG 
[20{ . the SF domains at incommensurate filling grow as 
the energy scale is lowered but there is no long-range su- 
perfiuidity because of one dimensional localization due to 
the disorder. This mechanism is very particular to one 
dimension, however, and does not directly carry over to 
higher dimensions. Whether or not commensurate fill- 
ings are special in 2D has to be tested with numerical 
simulations at non-integer filling. 

We note that the same type of anomalous exponential 
behavior ([2]) was previously observed in the magnetic sus- 
ceptibility of a quantum spin system with a glass phase 
similar to the MG [l^ . In that case a = 1 /2 was found in 
all cases studied, in contrast to the continuously evolving 
a we have found here in the BH model. One reason for 
the different a could be that z = 1 in the spin system. 
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